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Abstract. We study spectral properties of Dirichlet Laplacian on the conical layer 
of the opening angle n — 29 and thickness equal to tt. We demonstrate that below 
the continuum threshold which is equal to one there is an infinite sequence of isolated 
eigenvalues and analyze properties of these geometrically induced bound states. By 
numerical computation we find examples of the eigenfunctions. 
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1. Introduction 

Geometrically induced bound states in infinitely stretched regions with a hard-wall (or 
Dirichlet) boundary attracted a lot attention in recent years, in particular, because 
they represent a purely quantum effect which cannot be approached from the usual 
semiclassical point of view. The attention was first paid to tubes where the effective 



attraction is caused by bending |ES89[ IGJ921 IDE95j , later the analogous question was 
addressed with one dimension added, in curved layers [DEKOll ICEK04] . This question 
is more complicated and to the date we do not have a full analogue of the result saying 
that any local bending of an asymptotically straight tube exhibits such bound states. 
One reason is that geometry of surfaces is more complicated than that of curves and 
global properties play a more important role. 

However, this is not the only difference between tubes and layers. In the case of 
the former we have various solvable examples. They often correspond to tubes with 
a non-smooth boundary being thus beyond the assumptions under which the general 
results are proved and illustrating the robustness of the effect; on the other hand they 
can be composed of simple building blocks which facilitates their solvability. In fact, 



some examples came already before the general theory |LLM86l IESS891 ISRW89] and 
others were found in parallel with it — see, e.g., |LCM99] for a review. 

This is in contrast with our state of knowledge concerning the layers. We have 
reasonably strong and general existence results |DEK01t ICEK04j , later even extended 
to layers in higher dimensions [LL07] . as well a particular "weak coupling" result jEKOlj . 
but to the date no solvable example is knowii||. The goal of the present work is to discuss 
such an example in the geometrical situation which is particularly simple — apart of a 
trivial scaling the model depends on a single parameter. 

The layer to consider is the region bordered by two conical surfaces shifted mutually 
along the cone axis by vr esc 6. Using a natural separation of variables the angular 
component can be separated and the problem can be in each partial wave reduced 
to analysis of an appropriate operator on a skewed semi- infinite strip. This will be 
demonstrated in the next section and used in Sec. 3 to derive some properties of the 
corresponding discrete spectrum. To illustrate these results and to provide further 
insights into this binding mechanism, we present in the following section a numerical 
analysis of the problem: we calculate the eigenvalue dependence on the cone opening 
angle and present examples of eigenfunctions. 

I By solvability we do not mean finding the eigenvalues and eigenfunctions analytically, of course, 
this is beyond reach even in the tubes case. What we have in mind is the possibility of bringing the 
analysis to the state when one can find these quantities numerically and study their dependence on the 
parameters of the model. 
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Figure 1. The problem geometry. The left part shows the setting in cylindrical 
coordinates as described in Sees. I2.1] and l2.2l the right part shows the cross section in 
rotated coordinates s, u used in Sec. 12.31 



2. The Hamiltonian 

2.1. Basic definitions 

As we have said the object of our study is the Hamihonian of a quantum particle confined 
to a hard- wall conical layer of the opening angle tc — 29 for a fixed 9 G (0, ^tt), i.e. 

:= {(r, 0, z)eR': (r, z) e ilg, (I) e [0, 2tt)} , (2.1) 

where 

:= {{r,z) eR^ : r > 0, < z - r sin < vr esc 9} . (2.2) 

Without loss of generality we may assume that the layer width is vr, the result for 
other widths being obtained by a trivial scaling. Similarly we neglect values of physical 
constants and identify the Hamiltonian with the Dirichlet Laplacian — A^* on the 
Hilbert space "H = //^(Eg, rdrd(f)dz). Using the cylindrical coordinates we can write the 
operator action explicitly, 



i7V = - 



1 d 



dip 



1 d^ip d^ijj 



(2.3) 



r dr \ dr J d(f)'^ dz'^ ' 

with the domain consisting of all functions ip G W^^^ITjo) such that ip{x) = for 
X = (r, 0, z) G dTig and Hip G L^. The operator is associated with the quadratic form 
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(2.4) 



defined on functions ip G Wlp^ijle) vanishing at the layer boundary and such that 
(^^['j/'] < oo. It is well known |RS78l Sec. XIII.15] that the operator H is self-adjoint 
and positive, cr(if) C [0,oo). 



2.2. Partial wave decomposition 

Using the cylindrical symmetry of the system we write state Hilbert space as 
H = L^iVle, rdrdz) ® L^{Si) = L'^i^e, rdrdz) , 



(2.5) 
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where L'^{Si) refers conventionally to functions on the unit circle with the orthonormal 
basis {^e^'^'^ : m G Z}. The operator decomposes accordingly as 



m 



(2.6) 



with the domain consisting of G W^oci^e) satisfying the requirement ilj{r,r sin6) = 
ip{r, r sin ^^ + 7r tan 6) = and such that Hmip G L^. Under these conditions the operators 
H^, m 7^ 0, are self-adjoint while for m = one has to demand additionally that 
^|^_o ~ 0- boundary values here are understood as the appropriate limits. The 

quadratic forms associated with the partial wave components of the Hamiltonian are 
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(r, z) rdrdz (2-7) 
) = ■?/'(r, r sin 6' + tt csc 6') = and 



being defined on all ip G Wl^liVto) satisfying ip{r, r sin 6 
such that QmlV"] < 

As usual we can pass to unitarily equivalent operators by means of the trans- 
formation U : L'^{^lg, rdrd^;) — )■ L^(i7e) defined by {U'tp){r,z) := r^^'^ip(r,z) = ip{r,z) 
and an analogous map //^(Sg, rdrd^d^) — > L^(Se). In particular, the operator 
Hi := UHtU-^ acts as 
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(2.8) 



on functions 'ip G Wl^liyte) such that H^ip G satisfying the Dirichlet boundary 
condition at the boundary of VLq] the associated quadratic form is 



dip 



dr 



dip 
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+ 



m 



fr, z) drdz 



(2.9) 



with the domain of all ip G W^^^ {Vte) satisfying the Dirichlet condition at dVte and such 
that QmfV'] < oo. 



2.3. Alternative expressions 

There are other ways how to express action of our Hamiltonian. For instance, one can 
use rotated coordinates in the radial cross section plane according to Fig. [H 



s = r cos 6 + z sin 6* , u = —r sin 6 + z cos 6 , 
in which we have 

Vie ■= {{s,u) : s > 0, < m < min(7r, s cot 9)} . 
The partial wave operators act now as 

Am^ - 1 



ds^ du^ A[s cosO — usinOY 



(2.10) 
(2.11) 

(2.12) 
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the domain being the same as above; the associated quadratic form takes in these 
coordinates the shape 



ds 



du 



+ 



4m^ - 1 



4(s cos 6' — u sinOy 



One can also use skew coordinates, 

y = s — utanO , v = u, 
in which the expression f l2.13p acquires the form 



[s,u)dsdu. (2.13) 
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(2.15) 



3. Properties of the spectrum 

Since the conical layer is built over a surface with a cylindrical end its basic spectral 
properties follow from the theory developed on |DEK01l ICEK04j . 

Theorem 3.1. For any fixed 6 G (0, |7r) we have acss{H^) = [l,oo) while the discrete 
spectrum of the operator is non-empty with dim (Jdisc(-f^^) = oo. 

Proof. For most part of the argument we can refer to the papers quoted above. In 
particular, we proved there that inf cress(-f^^) = 1- In order to check that any k"^ + 1 
belongs to the spectrum one has to construct a suitable Weyl sequence. Consider three 
dimensional Cartesian coordinates s, t, u and take the family of functions 



9nis,t,u) :-- 
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f{ns, nt) cos 
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where \u\ < and / G C^(]R^) has the support in (—1,1)^. By construction 
||(-A - A;2 _ i)g^\\ ^ as n — )■ oo, and since the vanishing quantity is invariant 
w.r.t. Euclidean transformations we can find a family of functions cjn supported by 
rectangular blocks of the size 2r;, x 2n x (1 — ;|), each of them contained in Eg being 
naturally further and further from the layer tip as n grows, which proves the claim. 

The existence of the discrete spectrum was proven in [DEKOH ICEK04j . While it 
is not stated explicitly there that the number of the eigenvalues in (0, 1) is infinite, it 
is clear from the proof in which one constructs a sequence of trial functions supported 
by mutually disjoint ring-shaped domains of the conical layer, and thus orthogonal to 
each other; for large enough indices they give the value of the energy functional below 
the essential spectrum threshold. □ 
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After demonstrating the existence of the discrete spectrum let us look at its 
properties in more detail. First we note that all the geometrically induced bound states 
in the layer belong to the s-wave sector. 

Proposition 3.1. crdisc(-f^m) = ^ holds if m ^ 0. 

Proof. The claim is equivalent to crdisc(-f^m) — ^- For any higher partial wave, m 7^ 0, 
we have > /^^ |V^/'p(r, z) drdz which means that > —A^^ in the form sense 

where the rhs acts according to jRS78l Sec. XIII. 15] on functions ip ^ ^lof (^f) such 
that iplx) = for X = (r, z) G dQ^ and —Aip G L^. The result then follows by minimax 



principle because cr(— 



'-A 



D 



[1,00). 



□ 



This does not mean, of course, that the cylindrically asymmetric part of the problem 
is trivial — it can exhibit, for instance, an interesting scattering behaviour. Since we 
are dealing with the discrete spectrum here, however, we consider from now on only 
functions independent of the angular variable and drop the index m in the quadratic 
form expressions. By Theorem 13.11 the operator has a sequence of eigenvalues. 



(3.16) 



labelled in the ascending order; the corresponding eigenfunctions will be denoted as 
^0) g L^(S5i), j = 1, 2, . . . ,; with the symmetry in mind we use the same symbols for 
their radial cuts ip^^^ G L'^{VIq^ rdrdz) and ip^^^ for their counterparts in LF'iyLo). 

The main question is about the dependence of the spectrum on the cone opening 
angle. Let us start with demonstrating that the eigenvalues f l3.16p are smooth functions 
oi 9. To this aim we compare the cross sections (12. lip for two different values of 
the angle, 9 and 9' . They can be mapped onto each other by a longitudinal scaling, 
s cot 9 = s' cot 9'; then we can express the form ( 12.13^ referring to the angle 9' as 



tan^ 
tan 6" 



sin 29 



dip 
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[s, u) dsdu . 



(3.17) 



The integration is thus taken over the same domain and the angle dependence is moved 
to the coefficients of the operator. 

Proposition 3.2. The functions Aj(-), j = 1, 2, . . . , are on (0, ^tt). 

Proof. To find the derivative of Aj(-) we have Feynman-Hellmann theorem differentiate 
fIXTTj) with respect to 9' at 9' = 9. It yields 



sin 2^ 



ds 
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sin 2^ 



j) 



du 
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|^(^')|2cot2^ 
2(scos9 — usm9)'^ 



fs, u) dsdu . 
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Since ip'^^^ G Wy^^{Q,g) — in fact, it belongs to W^^^ — the first two contributions to 
tfie integral are finite. Combining this with the fact that Qm[ip^^^] < oo we infer that 
also the last contribution to the integral is finite, so the derivative exists. Using then 
the integrand continuity with respect to 6 and the dominated convergence theorem, we 
arrive at the result. □ 

Note that the terms appearing in the last expression have different signs so it is 
not apriori clear whether the eigenvalue behave monotonously as the angle changes. 
Nevertheless, the geometrically induced binding becomes stronger as 6 increases, in 
particular, we have the following result. 

Proposition 3.3. Set Aq := 7i"~^Joi ~ 0.58596 where jo,i is the first zero of Bessel 
function Jq. To any N ^ N and A G (Aq, 1) there is a 6q such that has at least N 
eigenvalues below A for any 9 G {9o, ^n). 

Proof. By Dirichlet bracketing |RS78t Sec. XIII.15] it is sufficient to find a region S C S^i 
such that — has this property. Consider a cylinder of radius R < n and length L. 
The eigenvalues of the corresponding Dirichlet Laplacian are easily found; it yields the 
condition 



We embed the cylinder axially into with one lid touching the tip of the inner surface; 
when 9 becomes close enough to we can chose R close to vr and L large so that the 
above condition is satisfied. □ 



4. Numerical results 



In order to get a more detailed information and a better insight let us study next the 
spectrum numerically. As a consequence of Proposition 13. H we may limit ourselves to 
two-dimensional geometry, i.e. to consider the spectrum of (12.61) for m = 0. Recall 
that the corresponding eigenfunctions are independent of the angular variable so we 
can choose = arriving at the two-dimensional problem depicted in Fig. [H Note also 
that the alternative formulation fl2.12p leads to the same results. 

The first difficulty we meet is the absence of reliable estimates concerning the span 
of eigenfunctions. The knowledge of the region where the eigenfunctions are negligible 
is important. We have to carry out integration over an infinite region and we need 
to cut the integration somewhere, unless we transform the problem to a finite region. 
This, however, brings problems of other kinds ruling out, in particular, use of the mode- 
matching approach. 

From general considerations we know that the smaller the gap between the 
eigenvalue and the continuum threshold, of course, the farther is the wavefunction 
stretched. Also, for small opening angles eigenfunctions spread farther. This rule of 
thumb may be of help as a rough estimate of integration region. 



Spectrum of Dirichlet Laplacian in a conical layer 



8 




Figure 2. The left figure depicts the reduced 2D region on which the solution is 
sought. The right one shows the dependence of eigenvalues on the opening angle; the 
first seven of them are shown (colour online). The essential spectrum threshold E — 1 
is the accumulation point of all of them and they approach it from below as — > k'^^- 
The thick dot marks the value Ao indicated in Proposition 13.31 the picture suggests 
that the eigenvalues tend to it monotonously as 9 ^ 0+. 



4.1. The method used 

We must admit that we have not found a method which would be efficient enough 
and intuitive at the same time. To summarize our experience we prefer finite element 
methods. In view of the absence of knowledge of a practical span of configuration space 
mentioned above even this approach requires separate treating of respective levels. In 
order to minimize the number of elements we took further advantage of symmetry and 
limited the problem to the region shown in the left part of Fig. [2] with the appropriate 
boundary conditions. The sought eigenvalues of the conical layer are then eigenvalues 
of the operator (12. 6p with m = on this region. 

We used 2nd order elements. The total number of degrees of freedom ranged from 
50000 to about 300000 depending on eigenfunction extent. All runs for small 6 (most 
demanding) took less than one minute each on an Intel Core 2 (1.86 GHz). 

We have also tried to employ the Galerkin method. The idea was to choose a set 
of test functions satisfying the boundary conditions. This, however, leaves still leeway 
in representing the subspace to which the problem is projected, i.e. how to choose the 
particular form of the test functions. The choices we attempted led to substantially 
more time consuming calculation than FEM, so fitting test functions are to be found 
yet. 
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4-2. The results 

The first question concerns the dependence of eigenvalues on the opening angle 
illustrated on Fig. [2] in the interval between 1 and 15 degrees. It agrees well with 
the general picture: the planar layer, 26 = it, does not support any eigenvalue below the 
threshold, while for any 6 < 71/2 there are by Theorem 13. II infinitely many such isolated 
eigenvalues, albeit all are weakly bound for sufficiently open conical layers. As the 
opening angle becomes more acute the low-lying states become gradually more strongly 
bound. The numerical results show no degeneracies. 

As we move away from the ground state, structure of the eigenf unctions acquires 
more and more a one-dimensional character. This is demonstrated in Fig [3] when the 
seven lowest eigenfunctions for the angle 6 = 2.5° are calculated. As expected, the 
ground-state eigenfunction can be chosen positive. On the other hand, the excited 
states display a number of nodal lines (nodal surfaces of the full problem) that is in 
a one-to-one correspondence with the eigenvalue index of the sequence {\i{6)}, in an 
analogy with the classical oscillation theorem. Notice also that the node distances are 
increasing with the distance of the cone tip. This behaviour has to be associated with 
the fact that the effective-curvature induced potential in (12.121) is proportional to 
having thus the same scaling behaviour as the kinetic part of the Hamiltonian. 

Another observation worth commenting is shown on Fig. HI where we have plotted 
(a side view of the) squared eigenfunctions l^/^^p with i = 1, . . . ,7 for the same opening 
angle 6 as above. As the index i grows, dominant regions lie around the inner tip and 
at the "far end" of eigenfunction extent, more exactly, between the last two nodal lines. 
Consequently, location of the regions with the largest probability of finding the particle 
in the conical layer depends on the eigenvalue number in this way. 

5. Concluding remarks 

Our discussion has not exhausted by far all open questions concerning spectra of conical 
layers. We have also mentioned the eigenvalue monotonicity which is expected to hold 
and supported by the above numerical result. A related question concerns the simplicity 
of the spectrum. It is again seen in the numerical example and its validity is strongly 
supported by the quasi-one- dimensional character of the higher eigenfunctions, but these 
considerations cannot replace a rigorous proof. 

Other open questions are associated with the behaviour for extreme values of 6 
which we left out in the previous section because these cases are numerically demanding. 
For small values of 6 the density of eigenvalues above Aq ~ 0.58596 should grow in 
accordance with Proposition 13.31 What we see nevertheless is that already for 6 = 2.5° 
the excited states have nodal lines (or nodal surfaces of the full problem) inside the cone 
"cap" ; one can ask about critical values of 6 for which one of those touches the inner tip. 
On the other hand, one can ask about the spectral asymptotics as ^ — ?■ ^tt— . In contrast 
to |EK01] the geometric perturbation is not compactly supported for conical layers and 
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Figure 3. The contour plot of the first seven eigenfunctions (colour online) for 
6 — 2.5°. In view of the axial symmetry, we show only one half of ip^r^z). The 
the r— and z— scales are different, the vertical axis being scaled by factor five: the 
transversal width is tt and z ranges from to 225. Nodal lines are indicated in black 
colour; higher eigenfuctions have clearly a quasi-onedimensional character. 
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Figure 4. Side view {i.e., with zero elevation angle) of lipl"^ for the same situation as 
in Fig. 3. All the eigenfunctions are normalized to the same constant, so that relative 
height ratios reflect the true situation (not distorted by different normalization of 
individual functions). There are two areas with higher probability, one in the vicinity 
of the inner tip and the other one at the "far end" of the function extent. The latter 
moves away from the vertex as the eigenvalue number increases. 
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the effective potential is on the borderhne between short- and long-range cases, so it is 
not a priori clear what the asymptotic expansion could be. 

Of interest is also the continuous spectrum of conical layers, both from the point of 
view of scattering — existence of wave operators and their asymptotic completeness — 
as well as of resonances. Being bound by the deadline of this issue, we leave all these 
questions to a later investigation. 
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